STATISTICAL PROPERTIES OF LYMAN BREAK GALAXIES AT z~4 
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Q H . Abstract. We study the luminosity function and the correlation function of about 1200 z~4 

Lyman break galaxies (LBGs) with i' < 26 that are photometrically selected from deep BRi' 
imaging data of a 618 arcmin 2 area in the Subaru/XMM-Newton Deep Field taken with Subaru 
Prime Focus Camera. The contamination and completeness of our LBG sample are evaluated, 
on the basis of the Hubble Deep Field-North (HDF-N) objects, to be 17% and 45%, respectively. 
We derive the UV (rest 1700A) luminosity functions (LFs) and find a large population of UV- 
luminous galaxies at z ~ 4. The LFs of the red and blue subsamples imply that the bright 
CN \ LBGs are redder in the UV continuum than the average color of the LBGs. Then we calculate 

the correlation function over 6 = 2" — 1000" and find that it is fitted fairly well by a power law, 



CO 

C3 



oj(6) = A^O , with A u = 0.71 ± 0.26. We estimate the correlation length ro (in comoving 



units) of the two-point spatial correlation function £(r) = (r/ro) -1 ' 8 to be tq = 2.7lg] 6 h~ 1 Mpc 
(Q m = 0.3 and Qa = 0.7). The correlation function shows an excess of to(0) on small scales 
(6 ^ 5"), departing from the power-law fit at > 3a significance level. Interpreting this as being 
due to galaxy mergers, we evaluate the fraction of galaxies undergoing mergers to be 3.0 ±0.9%, 
Q-f which is significantly smaller than those of galaxies at intermediate redshifts. 
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1 Introduction 



In the past several years, deep field surveys have made remarkable breakthroughs in studies on high 
redshift galaxies. The most successful survey would be the Hubble Deep Fields which are deep enough 
to study high-z galaxies up to z ~ 6. Many efforts on spectroscopic and photometric follow-up 
observations reveal hundreds of high-z galaxies and their properties (e.g.||). However the surveyed 
areas of HDFs are as small as a few arcmin 2 each, so they may be affected by field-to-field variations. 
Steidel and his collaborators have pioneered in statistical studies of high-z galaxies based on ground- 
based wide survey data |2(J. They isolate high-z galaxies in a two-color plane using the UV continuum 
features, and galaxies selected in this way are called Lyman break galaxies (LBGs). They studied the 



number density [21], the spatial distribution Q, nebular emission lines [|18|, and the stellar population 
[p!q]. Most of their studies are based on z ~ 3 LBG samples which were selected using U n GR colors. 
In this contribution, we extend the study to z ~ 4 LBGs based on our BRi' data taken with a 
newly installed Subaru Prime Focus Camera (Suprime-Cam: [i~2[] ) , which is a unique wide-field imager 
mounted on the 8m Subaru Telescope. Throughout this contribution, magnitudes are in the AB 
system, and all calculations assume a A-dominated spatially flat cosmology, (Q m ,Cl\) = (0.3,0.7). 



2 Observations and Data Reduction 

Deep and wide-field B-,V-,R-, and i'-band imaging data of a central 30' x 24' area in the Subaru/XMM- 
Newton Deep Survey Field (2 ft 18 m 00 s ,-5 o 12'00"[J2000]) were taken with Suprime-Cam during the 
commissioning observing runs on 2000 November 24-27. The present work is based on the B,R, and 
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Figure 1: Left panel: B — R vs. R — i' color diagram displaying the colors of model galaxies and stars. 
Typical spectra of elliptical, Sbc, Scd and irregular galaxies Q are redshifted from z = to z = 3, 
which are shown by four dashed lines. Each line has filled circles at z=0,l,2 and 3. We use a typical 
spectrum of galaxies at z ~ 3 to compute the z ~ 4 galaxy track which is shown by the solid line [13] . 
Star marks are 175 Galactic stars given by 0]. The arrow expresses the direction and magnitude of 
dust extinction on the B — R vs. R — i' plane for E(B — V) = 0.15, when Calzetti's extinction law 
H is applied. The box surrounding the upper left region is the selection criteria of our z ~ 4 galaxies 
(equation |l|). No obvious contaminant is found within the box. Right panel: B — R vs. R — i' color 
diagram displaying the colors of 5 a detected objects in our data. The box for our selection criteria is 
also shown. 

i' data. The individual CCD data were reduced and combined using IRAF and the mosaic-CCD data 
reduction software developed by us p3|] . The final images cover a contiguous 618 arcmin 2 area with a 
point-spread function FWHM of 0".98. The net exposure times of the final images are 177, 58, and 45 
minutes for B, R, and i' , respectively. The limiting magnitudes are B = 27.6, R = 26.5, and i! = 26.2 
for a 3<T detection in a 2" diameter aperture. The i'-band frame is chosen to detect objects, and we 
limit the object catalog to i! <; 26, in order to provide a reasonable level of photometric completeness. 



3 Selection of z~4 Lyman Break Galaxies 

Our catalog contains 42,557 objects with i' < 26.0 in total. On the basis of expectations from 
GISSEL96 Pj population synthesis models, we define the photometric selection criteria for galaxies at 
z ~ 4 (||13| for model parameters) as 

B-R> 2.1, R-i' < 0.5, B — R> 5.4(22 - »') + 0.9 (1) 

We estimate the redshift distribution, N(z), of galaxies satisfying equation ([!]) from the HDF-N objects 
for which magnitudes, colors, and photometric redshifts are given in j5|. We find that the criteria select 
z = 3.8 ± 0.5 galaxies. There are a total of 1192 objects that meet the criteria. Figure [l] shows the 
B — R vs. R — i' color diagram for model galaxies and Gunn & Stryker's stars M (left panel), and for 
5<r detected objects in our data (right panel). The left panel demonstrates that z ~ 4 galaxies are well 
isolated from interlopers, i.e., low- 2 galaxies and Galactic stars. We have estimated the contamination 
and completeness of our LBG sample by Monte Carlo simulations, generating artificial objects which 
mimic the HDF-N galaxies and distributing them on our original images. The contamination is defined, 
for the detected simulated objects, as the ratio of low-redshift (z < 3.3) objects meeting equation (|l|) 
to all the objects satisfying equation ([!]). The completeness is defined as the ratio of z — 4 simulated 
objects passing our detection threshold and satisfying equation ([[]) to all (detected + undetected) 
z ~ 4 simulated objects. We find from the simulations that the completeness and the contamination 
are 45% and 17%, respectively. 
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Figure 2: (a) Luminosity function (LF) of z ~ 4 LBGs, together with those for z = and 3 galaxies. 
The filled circles indicate z ~ 4 LBGs derived from our data, and the squares are for the z ~ 4 
galaxies from the HDF photo- z catalog ||. The solid line shows the best fitted Schechter function to 
the z ~ 4 data. The dashed line is the LF of UV selected galaxies at z = (whose magnitudes are 
in Mab (2000A) |2|].) The dotted line denotes the LF of z ~ 3 LBGs A large population of 

UV-luminous (Mab —20) galaxies is seen at z ~ 3 and 4. (b) LFs of the red and blue subsamples. 
The filled circles and filled squares show the red sample, while the open circles and open squares are for 
the blue sample. The solid line is the best fitted Schechter function obtained from the whole sample, 
which is shown as the solid line in (a). See text for details. 



4 Luminosity Function 

The UV luminosity function (LF) of z ~ 4 LBGs is derived from our 1192 objects with i' < 26. We 
calculate the surveyed effective volume which is the integrated volume from z = 3.3 to z = 4.3 by 
taking into account the completeness and contamination of the sample selection, and N(z) (§||). The 



LF is shown in figure |2|(a), together with those at z = [22] and z ~ 3 [21]. The LF of our z ~ 4 
LBG sample is consistent with the one derived by [21]. We fit the Schechter function to the LF, and 



find M^ B (1700A) = -20.41, (p* = 1.9 x 10" 3 with a fixed slope a = -1.6. It is found from figure 
H(a), that all the LFs from z=0 to 4 seem to prefer a steep slope, a ~ —1.6. The M*s at z ~3 and 
4 are about 2 magnitude brighter than that at z=0. The number of bright (Mab (1700^4) <J —20) 
galaxies are much larger at high redshifts than in the local universe, while that of faint galaxies are 
comparable. This implies that a large population of UV-luminous galaxies exists only at high redshifts. 
This population reflects the high cosmic starformation rates at high redshifts, which were derived by 
[ pX| 1 and |2l|| . There are two extreme explanations for this population. One is that the star-formation 
rate of z ~ 3 — 4 galaxies is intrinsically higher than that for z = galaxies, and the other is that the 
z ~ 3 — 4 galaxies are very dust-poor (or partially dust free) and thus the M* values are apparently 
brighter than that for z = galaxies which suffer from non-neglibible dust extinction. A detailed 



discussion will be given in [14]. 



We divide our whole LBG sample into two subsamples, blue and red, by the R — i' color. The 
blue (red) sample is composed of objects bluer (redder) than the median color of the whole sample, 
(R — i') = 0.12. The LFs of red and blue samples are calculated in the same manner as for the whole 
sample. Figure ^(b) shows the LFs of red and blue samples. These LFs are not significantly different 
at the faint end (Mab > —20.3), however, the LF of red sample seems to be higher than that of the 
blue sample at the bright end (Mab < —20.3). This may imply that the bright LBGs are dominated 
by galaxies red in the R — i' color. The R — i color reflects the shape of UV continuum between 
1380 — 1620A for z = 3.8 LBGs Q and the difference in the color is thought to be caused by various 

1 1530A-1810A for z = 3.3 LBGs, 1250A-1470A for z = 4.3 LBGs 
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Figure 3: (a) Sky distribution of the z ~ 4 LBGs. Different symbols correspond to different magnitude 
bins defined in the panel. Masked regions, to avoid the effects of bright stars, are shown as dotted 
circles, (b) Angular correlation function for the z ~ 4 LBGs. The solid line shows the best-fit power 
law with iv (8) = A u 
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reasons, for example, differences in dust extinction, stellar population, and redshift. If the difference 
in the color is mainly originated from the dust extinction, the boundary of the color, R — i = 0.12, 
is estimated to be E(B — V) ~ 0.18, following [11]. If this is the case, the brighter LBGs may be 
dominated by dusty LBGs with E{B — V) > 0.18. Note that in deriving the red and blue LFs, we 
assume that the completeness and the contamination of these subsamples are the same as for the 
whole sample. We should point out the possibility that the difference in red and blue LFs comes from 
a systematic difference in the completeness and contamination factors of the blue and red samples. 



5 Angular Correlation Function 

Figure |3](a) shows the sky distribtuion of z ~ 4 LBGs in our sample. In this figure, we find a somewhat 
inhomogeneous distribtution of LBGs, especially for bright (i! < 24.5) ones. We derive the angular 
two-point correlation function u(8) using the estimator defined by || (see Jl3[] ). The resulting angular 
correlation function for the sample is shown in figure ||(b) 0. 

We fit a power law, (v(8) = A U 8~P , to the data points, and find the slope j3 = 0.6^q]|, which is 
consistent with those for nearby galaxies, = 0.7 — 0.8, though the errors in our estimate are large. 
The best fit value for j3 = 0.8 is A w = 0.71 ± 0.26 in units of arcsec^. 

The two-point angular correlation function is related to the spatial correlation function £(r) 



(r/ro) 7 by an integral equation, the Limber transformation [17] . We apply the Limber transformation 



to the best fit value of with (3 = 0.8, using the redshift distribution N(z) of the LBGs (pH. The 



obtained correlation length is tq = 2.1~^^ i h~ 1 Mpc. This is twice as large as the correlation length 
of the dark matter at z ~ 4 predicted from an analytic model by flq| . The difference is thought 
to be due to the biasing of galaxy distribution, which is characterized by the biasing factor, b, as 
igai{r) = b 2 '^matter ■( r )- We obtain a linear bias b ~ 2.6, which is similar to that found by for bright 
z ~ 3 LBGs. On the other hand, the semianalytic models predict much stronger clustering for galaxies 
(eg. |[|], ||). A further discussion requires a detailed comparison of galaxy properties and selection 
criteria, between galaxies in the models and LBGs in our sample. 

In figure |||(b), we find an excess ^ of tv(8) at small scales (8 <; 5") at > 3a significance levels, 

2 Error bars do not include sample variance caused from field-to-field variations. The effect of the variations in our 
sample is propably modest since our sample probes a large comoving volume, 36 x 36 x 520 = 6.7 x 10 5 h _1 Mpc 3 . 

3 With a visual inspection, we find that all the galaxies contributing to the excess are real; no artificial objects, such 
as halos of bright stars, are found to be included in them. 



relative to the best-fit power law. This can be caused by various effects; for example, galaxy-galaxy 
mergers, two galaxies within a common dark matter halo, and/or field variance. If we interpret this as 
galaxy-galaxy mergers of those i' < 26 LBGs, which are classified as bright LBGs with M(1700A) <^ 
M*(1700A) (§||), we can calculate, following jH]], the fraction of galaxies undergoing mergers as 
from the observed and expected numbers of galaxy pairs with a separation of l'/5 < 9 < 4" (the 
projected physical separation corresponds to 7.5 to 20 7i kpc at z = 3.8). The resulting value is 
fpair = 3.0 ± 0.9%. This value is significantly smaller than those for intermediate-redshift {z ~ 0.3) 
galaxies whose f pa i r at <^ 20 /i _1 kpc is 5 — 15% (e.g. £L9|). 
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